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ABSTRACT 

HATS-8b is a low density transiting super-Neptune discovered as part of the HATSouth project. The 
planet orbits its solar-like G dwarf host (V=14.03 =b 0.10, T e ^=5679 =b 50K) with a period of 3.5839d. 
HATS-8b is the third lowest mass transiting exoplanet to be discovered from a wide-held ground 
based search, and with a mass of 0.138 d= 0.019 Mj it is approximately half-way between the masses 
of Neptune and Saturn. However HATS-8b has a radius of 0.873^0.075 resulting in a bulk density 
of just 0.259 =b 0.091 gem -3 . The metallicity of the host star is super-Solar ([Fe/H] =0.210 ± 0.080), 
arguing against the idea that low density exoplanets form from metal-poor environments. The low 
density and large radius of HATS-8b results in an atmospheric scale height of almost 1000 km, and 
in addition to this there is an excellent reference star of near equal magnitude at just 19 /r separation 
on the sky. These factors make HATS-8b an exciting target for future atmospheric characterization 
studies, particularly for long-slit transmission spectroscopy. 

Subject headings: planetary systems — stars: individual (HATS-8) — techniques: spectroscopic, 
photometric 
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This paper includes data gathered with the 6.5 m Magellan Tele¬ 
scopes located as Las Campanas Observatory, Chile. Based in 
part on observations made with the MPG 2.2 m Telescope and 
the ESO 3.6 m Telescope at the ESO Observatory in La Silla. 
This paper uses observations obtained with facilities of the Las 
Cumbres Observatory Global Telescope. 


While the Solar System planets display a rich diver¬ 
sity of physical properties, the discovery of exoplanets 
over the last two decades have revealed an even wider 
range of systems. The mass regimes of exoplanets is 
an example of this expanding d i versit y. The discovery 
of “super-Earths” (|Rivera et all I20051) . with masses of 
2 < M® < 10, have shown us a type of planet unlike any¬ 
thing in our Solar System. It now appears that planets 
in this mass range a re relatively common in the Galaxy 
(|Fressin et al.l l2013h . although for many planets discov¬ 
ered by Kepler the masses can only be inferred from radii 
in a statistical sense (using mass-radius relationships) 
due to the difficulties of measuring radial velocities or 
transit timing variations for small planets orbiting faint 
host stars. 

Another mass-class of planets that we only find outside 
our Solar System are planets more massive than Neptune 
at 0.05Mj, but smaller than Saturn at 0.3Mj. These low 
mass gas planets, or “super-Neptunes” are at the very 
transition between planets with H/He dominated com¬ 
positions and those without H/He dominating the bulk 
composition. In this paper we present the discovery of 
HATS-8b, a transiting super-Neptune found as part of 
the HA TSou th survey for southern transiting exoplanets 
IBakos et all (f20l! . HATS-8b has a mass lying almost 
half-way between Neptune and Saturn. One of the great 
advantages of ground-based transit surveys that target 
bright stars is that we can measure both the mass and 
radius of discovered exoplanets, and therefore we are able 
to determine the density of the exoplanet. This is par¬ 
ticularly important at these mass-ranges, as we cannot 
infer the radius from the mass alone or vice versa. 

In Section [2] of this paper we outline the photomet¬ 
ric and spectroscopic observations that led to the detec¬ 
tion of HATS-8b. In Section [3] we detail the methods we 
used to determine the physical parameters of the planet 
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and star as well as ruling out non-planetary interpreta¬ 
tions of the data. Finally, in Section [4] we put HATS- 
8b into context with other exoplanets discovered in this ra 
mass/density range, and discuss the fact that this exo- e 
planet orbits a super-solar metallicity star. We also dis¬ 
cuss the possible follow-up opportunities for this system 
from the ground and from space. 

2. OBSERVATIONS 
2.1. Photometric detection 

HATS-8 was in tensivel y monitored as part of the HAT- 
South survey dBakos et al.l l2013h . Although some data 
was acquired as early as September 2009, the bulk of 
the observations, over 10,000 images, were taken between 
March 2011 and August 2011 (see Tabled]). Details of the 
HATS o uth imaging system can be found in iBakos et alJ 
(1201311 ; iPenev et al.l (1201311 , while here we provide a sum¬ 
mary of the critical features. HATSouth employs Taka- 
hashi astrographs (//2.8, 18cm apertures) imaged on to 
Apogee U16M 4K x 4K cameras. Imaging is performed 
in the Sloan r band with exposure times of 240s and a 
mean cadence of approximately 300s. Images are col¬ 
lected from all three sites in the global network (see Ta¬ 
ble \T\ which also gives a breakdown of the number of 
images taken at each site), and data is reduced and light 
curves produced v ia apertu re photometry in the man¬ 
ner detailed in iPenev et al.l 1 2013 1. Light curves are de¬ 
trended u sin g Ext ernal Parameter Decorrelation (EPD; 

IBakos et al.l 201 0h a nd the Trend Filtering Algorithm 
(TFA; iKovdcs et 511200511 1. The detrended light curve 
for HATS-8 is set out in Table [2] We search light curves 
for periodic transit-like ev ents u sing the B ox-fi tting Least 
Squares algorithm (BLS iKovacs et ahl 12002( 1. For the 
case of HATS-8, we find a transit signal at a period 
of P=3.5839 days with a depth of 8.3 rnrnag (see Fig¬ 
ure [T]). No observable out-of-transit variation or sec¬ 
ondary eclipse is apparent. We therefore began recon¬ 
naissance spectroscopic observations of this candidate as 
is detailed in Section l2~2l 

2.2. Reconnaissance Spectroscopy 

Although a periodic dip in the light curve of a star 
may be due to a transiting exoplanet, there are many 
other astrophysical sources of such a signal, which can 
be termed “false positives”. Blended eclipsing binaries 
and grazing eclipsing binaries are common causes of false 
positives. Other candidates can be “false alarms”; that is 
the detected event is non-astrophysical in nature. These 
can be caused by artifacts on the CCD or period sys- 
tematics introduced by the telescope or mount systems. 

The combination of all of these false candidates, cou¬ 
pled with a desire to find shallow transits (i.e. smaller 
radius planets), means that the majority of candidates 
produced by ground-based transit surveys such as the 
HATSouth survey are not genuine transiting exoplanets. 

It is therefore crucial to undertake an efficient vetting 
program of transiting exoplanet candidates by way of 
reconnaissance spectroscopy before attempting the re¬ 
source intensive search for subtle radial velocity varia¬ 
tions consistent with an orbiting exoplanet. Such recon¬ 
naissance spectroscopy will bec ome even m ore cr itical 
as large surveys such as T ESS dRicker et al .11201 -H) and 
PLATO (|Rauer et a l. 2014) produce enormous numbers 
of candidates requiring such vetting. 
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Fig. 1.— The HATSouth discovery light curve for HATS-8b. In¬ 
dividual points show the unbinned instrumental r band light curve 
of HATS-8 folded with the period P = 3.5838933 days resulting 
from the global fit described in Section [3] The solid line shows the 
best-fit transit model (see Section The lower panel is a zoom-in 
on the transit feature, which has a depth of 8.3 mmag and duration 
of 3.031 hours. The dark filled points in the lower panel show the 
light curve binned in phase using a bin-size of 0.002. 

The candidate HATS-8b is a good example of low a 
signal-to-noise detection, with a transit depth in the dis¬ 
covery light curve of just 8.3 mmag, compared to an 
RMS precision of the light curve per observation of 12 
to 14 mmag level (see Table |T|. 

In the case of HATS-8, follow-up spectroscopy was per¬ 
formed in October 2012 using the echelle spectrograph on 
the 2.5 m Du Pont telescope at Las Campanas Observa¬ 
tory in Chile. We obtained two spectra using the 1" x 4" 
slit (R~40000) on the nights of October 25 and October 
26 2012, each with an exposure time of 1800 s. From 
these observations we calculated that HATS-8 was a G- 
dwarf with a low projected rotational velocity (usin*) 
and there was no sign of a secondary spectrum that could 
be indicative of a binary system. A further observation 
of HATS-8 was obtained with the FEROS spectrograph 
(iKaufer fc Pasquinlll998l ) on the MPG 2.2 nr telescope at 
the ESO Observatory in La Silla, Chile. This instrument, 
with slightly higher spectral resolution, confirmed the Du 
Pont finding that HATS-8 is a G-dwarf with low usini 
and no evidence of a composite spectrum. The three 
observations also showed no radial velocity variation to 
within the uncertainty of the measurements, indicating 
the system could not be a eclipsing binary system with a 
large radial velocity amplitude. Details of the three re¬ 
connaissance spectra taken for HATS-8 are set out in Ta¬ 
ble [3] Based on the reconnaissance spectroscopy, HATS- 
8 was deemed to be a sufficiently strong candidate to 
warrant undertaking precise photometric and radial ve¬ 
locity measurements as detailed in Sections 12.31 and 12.41 
respectively. 

2.3. Precise Photometric Follow-up 

Precise photometric observations of a transit of HATS- 
8b were carried out using the SiTe3 imaging camera on 
the SWOPE lrn telescope at Las Campanas Observatory 
(LCO) in Chile on 2013 May 29. Defocused imaging over 
the 14.8'x22.8' field was performed in i band with expo¬ 
sure times of 120s. The images were reduced and aper- 
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TABLE 1 

Summary of photometric observations 


Facility 

Date(s) 

Number of Images a 

Cadence (s) b 

Filter 

Precision (mmag) 

HS-1/G579 

2009 Sep-2011 Aug 

4484 

301 

r band 

13.9 

HS-3/G579 

2010 Mar-2011 Aug 

2607 

303 

r band 

13.3 

HS-5/G579 

2010 Sep-2011 Aug 

3303 

303 

r band 

12.2 

Swope/SiTe3 

2013 May 29 

107 

213 

i band 

2.3 

Swope/E2V 

2014 Jul 01 

116 

189 

i band 

1.8 

LCOGT lm/Sinistro 2014 Jui 09 

50 

347 

i band 

3.5 


a Excludes images which were rejected as significant outliers in the fitting procedure. 

k The mode time difference rounded to the nearest second between consecutive points in each light curve. Due to visibility, 
weather, pauses for focusing, etc., none of the light curves have perfectly uniform time sampling. 


ture photo metry ex t racte d in the standard manner as set 
out in iPenev et al.l (2013). Figure [2] shows the resulting 
light curve, which covers all but the ingress of the transit 
on that night, and allows for a very precise determina¬ 
tion of the period and phase of HATS-8b. Using this 
precise ephemeris, a second photometric observations of 
transits of HATS-8b was obtained with the 1 m SWOPE 
telescope on the night of 20f4 July 1. Again the moni¬ 
toring was carried out in i band, this time using the E2V 
imaging camera with exposure times of 150s. The egress 
of the transit of HATS-8b was detected with these obser¬ 
vations, as detailed in Figure [2] Finally, nine days later 
we observed another transit of HAT S-8b, t his time usin g 
the LCOGT lm telescope network dBrown et all 12013 ). 
This observation was carried out with the LCOGT lm 
telescope at CTIO in Chile using the Sinistro imaging 
camera in i band with a 300 s exposure time and the 
telescope slightly defocused. The observation covered a 
full transit. The transits are shown in Figure [2] 

All three transit events monitored with higher preci¬ 
sion photometry confirmed the transit signal detected in 
the HATSoutlr discovery light curves. The depth and 
shape of the transit are consistent with that expected 
from a transiting exoplanet, and we use these data in 
combination with the discovery photometry to determine 
the global parameters of the HATS-8b system in Sec¬ 
tion [3] 

2.4. High Precision Spectroscopy 

Radial velocity measurements play a critical role in 
the confirmation and characterization of transiting plan¬ 
ets, and are needed to determine the mass, and hence 
the bulk density, of the exoplanet. HATS-8 is a difficult 
prospect for radial velocity measurements, as it is rela¬ 
tively faint (V=14.03). For this reason we needed to use 
Keck/HIRES in order to have the necessary precision and 
signal-to-noise to measure the radial vel ocity variations . 

We observed HATS-8 with HIRES dVoet et al.l I1994T) 
on Keck-I in Hawaii between June and September 
2014. We used HIRES in its standard configuration 
for precise radial velocity measurements: a slit width 
of 0.86", A/AA ~ 55,000, and wavelength coverage of 
3800-8000 A. Exposure times were typically 1500 s and 
achieved a signal-to-noise ratio of 40 per pixel in the con¬ 
tinuum near 5500 A in the reduced HIRES spectra. An 
iodine gas absorption cell in the optical path is used to 
superimpose iodine absorption lines on the stellar spec- 
trum and provide an a ccurate wavelength calibration (see 
iMarcv fe Butleill 19921) . A synthetic iodine-free template 


spectrum is created by interpolating the ICoelhol (I2014D 
grid of stellar models for the values of T e f p, logg, [Fe/H] 
derived in Section l3Jl We found that using the syn¬ 
thetic template produced smaller residuals and measure¬ 
ment uncertainties when compared to using an observed 
iodine-free spectrum with moderate S/N as the template. 
In the case of the synthetic template, no deconvolution 
is needed to remove the effects of an asymmetric instru¬ 
ment al PSF, and it is completely noise-free. See iFultonl 
(1201511 for a more detailed description and performance 
tests of the synthetic template technique. We derived 
relative radial velocity measurement for HATS-8 using 
the method described in Butler et al. (1996), which ac¬ 
counts for variations of the spectrograph instrumental 
line profile. The radial velocity measurements are listed 
in Table 0] and plotted in Figure 01 along with the esti¬ 
mated uncertainties in these measurements and the bi¬ 
sector spans of the average spectral line. 

We computed spectral line bisector spans from the 
blu e orders of the Keck/HIRES observations following 
iTorres et al.l ([2007) and corrected these for contamina¬ 
tion fro m scattered moonlight following Ha rtman et al.1 
(l2011alf . These data are also represented graphically as a 
function of orbital phase in Figure 0J along with the best 
fit circular orbit. No systematic variation is seen in the 
bisector span measurements, which would have been a 
signature of a blended system. The overall scatter of the 
bisector spans is only 9.5 ms -1 . This is good evidence 
that HATS-8 is not a blended stellar eclipsing binary sys¬ 
tem. A more thorough rejection of blends is discussed is 
Section 13.21 

We measured the emission in the cores of the 
Ca II H&K lines and found a logi?[j K value of -5.179. 
This value suggests that this star has little chromospheric 
activity. The residuals of the radial velocity fit show 
an RMS of 5.1ms _1 , which is consistent with the for¬ 
mal uncertainties (see Figure [3]) and implies a very faint, 
chromospherically quiet star. 

3. ANALYSIS 

3.1. Stellar Properties of HATS-8 

The an alysi s of the reconnaissance spectra, described 
in Section !?^! indicated that the star HATS-8 is a slowly 
rotating G-dwarf. We derive more precise stellar param¬ 
eters by analyzing the iodine-free template spectra taken 
with Keck/HIRES on 2014 June 20 using the Stellar Pa¬ 
rameter Classification (SPC) m ethod as describe d in th e 
supplementary information to iBuchhave et al.1 (l2012f) . 
This method involves cross-correlating the observed spec- 
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TABLE 2 

Differential photometry of HATS-8 


BJD 

(2 400 000+) 

Mag a 

<T Mag Mag(orig) b 

Filter 

Instrument 

55774.34619 

0.00664 

0.00797 

r 

HS 

55699.08460 

-0.02530 

0.01660 

r 

HS 

55709.83639 

-0.00581 

0.02028 

r 

HS 

55770.76289 

-0.02414 

0.00927 

r 

HS 

55763.59560 

-0.00056 

0.00983 

r 

HS 

55788.68300 

-0.00119 

0.01462 

r 

HS 

55763.59611 

-0.01683 

0.01040 

r 

HS 

55691.91833 

-0.04584 

0.01747 

r 

HS 

55663.24726 

-0.01226 

0.00956 

r 

HS 

55727.75748 

0.02984 

0.02475 

r 

HS 


Note. — This table is available in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content. The data are 
also available on the HATSouth website at http://www.hatsouth.org 


a The out-of-transit level has been subtracted. For the HATSouth light curve (rows 
with “HS” in the Instrument column), these magnitudes have been detrended using 
the EPD and TFA procedures prior to fitting a transit model to the light curve. 
Primarily as a result of this detrending, but also due to blending from neighbors, the 
apparent HATSouth transit depth is somewhat shallower than that of the true depth 
in the Sloan r filter (the apparent depth is 85% that of the true depth). For the 
follow-up light curves (rows with an Instrument other than “HS”) these magnitudes 
have been detrended with the EPD procedure, carried out simultaneously with the 
transit fit (the transit shape is preserved in this process). 

k Raw magnitude values without application of the EPD procedure. This is only 
reported for the follow-up light curves. 


TABLE 3 

Summary of spectroscopic observations 


Telescope/Instrument 

Date Range 

Number of Observations 

Resolution 

Observing Mode 

du Pont 2.5m/Echclle 

2012 Oct 25-26 

2 

40000 

RECON 

MPG 2.2 m/FEROS 

2013 July 18 

1 

48000 

RECON 

Keck-I lOm/HIRES 

2014 June 20 

1 

55000 

I2-free template 

Keck-I lOm/HIRES 

2014 June 17 - 2014 Sep 10 

9 

55000 

I2/HPRV 


TABLE 4 

Relative radial velocities and bisector span measurements of 

HATS-8. 


BJD 

(2 456 000+) 

RV a 

(ms- 1 ) 

ffRV b 

(ms -1 ) 

BS 

(ms -1 ) 

CBS 

Phase 

Instrument 

826.94201 

-8.93 

6.15 

10.0 

16.4 

0.202 

Keck 

827.93901 

-2.39 

6.24 

-12.6 

13.5 

0.480 

Keck 

828.95001 

17.25 

7.19 

-7.0 

10.8 

0.762 

Keck 

830.95801 

-25.03 

6.34 

4.7 

5.5 

0.323 

Keck 

874.03501 

-19.55 

12.68 



0.342 

Keck 

882.81701 

13.43 

6.20 

-2.5 

11.0 

0.793 

Keck 

889.79501 

17.38 

5.65 

4.5 

5.9 

0.740 

Keck 

891.80901 

-21.07 

6.75 

8.6 

9.6 

0.302 

Keck 

911.74301 

14.49 

6.23 

-14.8 

7.4 

0.864 

Keck 


a The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted 
separately to the Keck/HIRES velocities in Section [3] has been subtracted, 
k Internal errors excluding the component of astrophysical/instrumental jitter 
considered in Section [3] 
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Fig. 2. — Left: Follow-up % band transit light curves of HATS-8 from the Swope lm and LCOGT lm telescopes in Chile. The dates 
and instruments used for each event are indicated. The light curves have been detrended using the EPD process. Light curves after the 
first are shifted for clarity. Our best fit is shown by the solid lines. Right: Residuals from the fits in the same order as the curves at left. 
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Phase with respect to T c 

Fig. 3. — Top panel: High-precision radial velocity measure¬ 
ments for HATS-8 from Keck/HIRES, together with our best-fit 
circular orbit model. Zero phase corresponds to the time of mid¬ 
transit. The center-of-mass velocity has been subtracted. Second 
panel: Velocity O — C residuals from the best-fit model. The error 
bars for each instrument include the jitter which is varied in the fit. 
Third panel: Bisector spans (BS), with the mean value subtracted. 
Note the different vertical scales of the panels. 

trum between 5050 and 5360 A with a grid of synthetic 
templates covering a wide range of T e g, vsin*, log g, and 
[Fe/H], 

Using the method described in iSozzetti et al.l (|2007t l 
and app lied in previous HATSouth discoveries (e.g. 
iPenevet alJ 12013 ). we initially derive the mean stellar 
density of HATS-8 via light curve fitting, which we de¬ 
termine to be 1.20±o 29 g cm 3 . To determine other stel¬ 
lar properties, we combine this mean stellar density with 
the Tpff ^ from SPC analysis, and the Yonsei-Yale (Y2; 
lYi et all 12001 1 stellar evolution models. This yields a 
logg*=4.40±0.06, which we then fix in a second itera¬ 
tion of SPC to derive the final parameters for HATS-8. 
The final parameters are listed in Table 0 and reveal 
HATS-8 to be G dwarf host with a T e fj=5679 ± 50 K, a 
log(7=4.386±0.07i, and a v sin *=2.00 ±0.50 km s _1 . We 
note that this stellar surface gravity (log g=4.386±0.07i) 
is consistent with the value measured from spectra alone 
(log5=4.48 ± O.fO). We estimate the age of this system 
using using the Y2 model isochrones to be 5.1 ± 1.7 Gyr; 
see FigurelU All of these stellar parameters make HATS- 
8 very similar to our own Sun. The primary differ¬ 
ence is that the metallicity of HATS-8 is super-Solar at 
[Fe/H]=+0.210± 0.080. 

3.2. Excluding Blending Scenarios 

One of the most difficult false positive scenarios to rule 
out is a background blended eclipsing binary, which is 
not obvious in either photometric or spect roscop ic follow¬ 
up. He re we follow the method set out inlHartma n et all 
(I20IH) in order to determine the if the data can be ex¬ 
plained by a blended stellar system rather than a tran- 



Fig. 4.— Comparison between the measured values of T e ff* 
and p* (from SPC applied to the HIRES spectra, and from our 
modeling of the light curves and radial velocity data, respectively), 
and the Y 2 model isochrones from lYi et all 1200 il l. The best-fit 
values (dark filled circle), and approximate lcr and 2a confidence 
ellipsoids are shown. The values from our initial SPC iteration are 
shown with the open triangle. The Y 2 isochrones are shown for 
ages of 0.2 Gyr, and 1.0 to 14.0 Gyr in 1 Gyr increments. 

siting exoplanet. 

We find that a model consisting of a single star with 
a transiting planet provides a lower x 2 fit to the pho¬ 
tometric data than any of the blended eclipsing binary 
models that we tested. The best-fit blended eclipsing bi¬ 
nary model can be rejected in favor of the star+planet 
model, based solely on the photometry, with greater than 
2cr confidence. We also simulated cross-correlation func¬ 
tions (CCFs), radial velocities and bisector span mea¬ 
surements for each blended eclipsing binary model con¬ 
sidered, and find that all of the models that come close 
to fitting the photometric data would have easily been 
detected as composite stellar systems based on the spec¬ 
troscopic observations (they have clearly double-peaked 
CCFs and/or radial velocities and bisector spans that 
vary by more than 1 km s _1 ). We therefore conclude that 
HATS-8 is a transiting planet system. 

We also considered the possibility that HATS-8 is a 
transiting planet system with an unresolved stellar com¬ 
panion next to the host star. While we cannot rule out 
very low mass stellar companions, we find that a compan¬ 
ion with M > 0.65 M e can be excluded with greater than 
3cr confidence. If HATS-8 has a 0.65 M 0 companion, the 
radius of HATS-8b would be ~ 6% larger than what we 
measure assuming no stellar companion; the mass would 
also be slightly larger. 

3.3. Global Modeling of Data 

We derive the parameters for the HATS-8b system us¬ 
in g a j oin t, Ma r kov-C hain Monte Carlo method detailed 
in|Bakos et al.| (12010( 1. and with modifications set out in 
lHartman et al.1 ( 2012 ). This method models the best fit 
parameters in a global way using the discovery photom¬ 
etry, follow-up photometry, and precise radial velocity 
measurements. The resulting planetary parameters are 
set out in detail in Table Ea We model the system for 
both a circular planet orbit and for an orbit where ec¬ 
centricity is allowed to vary. We calculate the Bayesian 
evidence for both these scenarios, and find that the cir¬ 
cular orbit is preferred. However we also list the 95% 
confidence upper limit, which is e< 0.376. 

Our global modeling reveals that HATS-8b is a 0.138± 
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0.019 Mj exoplanet, which is just 2.5 times the mass 
of Neptune or one-third the mass of Saturn. The ra¬ 
dius is 0.873lo.o75 w hich results in a bulk density 
of just p p =0.259 ± 0.091 gem -3 . Its orbital distance of 
0.04667 ± 0.00055 AU implies that HATS-8b would have 
a temperature of 1324±g| K (assuming zero albedo and 
a complete redistribution of heat). 


4. DISCUSSION 
4.1. Density 


HATS-8b is a very low density exoplanet, with a mean 
density of just 0.259 ± 0.091 gem -3 . Figure [5] shows 
HATS-8b in context with all known exoplanets with pre¬ 
cisely measured densities. We see from Figure [5] that 
HATS-8b sits in a transition region: between the low- 
mass, non-degenerate planets where bulk density de¬ 
creases with total mass and the high mass, partially de¬ 
generate gas giants where bulk density increases with to¬ 
tal mass. Discovering transiting super-Neptunes in this 
transition region is important, as it allows us to investi¬ 
gate the properties of pla nets that ha ve n ot undergone 
“run-away” gas-accretion (iMordasini et al.ll2009H , in con¬ 
trast to the well-studied population of hot-Jupiters. 

HATS-8b is close st in mass and densit y to the exo¬ 
planet HAT-P-18b (IHartman et al.lfeOllaf) . which has a 
similar density (0.25 g cm~ 3 ) bu t is slightly mo_re massive 
at M=0.197Mj. Using the iFortnev et all (1200711 models 
for gas giant planets, and assuming an age of the system 
to be approximately 4.5 Gyrs as indicated from the YY 
isochrones (see Figure |4}, HATS-8b would have a core¬ 
mass of just IOMq. As with HAT-P-18, the metallicity of 
HATS-8 is not low ([Fe/H]=0.210±0.080), and therefore 
the low-mass core that we infer is unlikely to be related to 
a lack of metals in the proto-planetary disk from which 
HATS-8b formed. This argues against the metallicity- 
ra dius relationship f or low-mass gas giants as discussed 
in lFaedi et alJ (1 2011 :). Purely in terms of mass, HATS-8b 
is similar to the exoplanet Kepl er-1 0 1b, which is slightly 
more massive at 0.16Mj dBonomo et al.l l2bl41. However 
Kepler-lOlb has a radius of just 0.515i?j and a resulting 
bulk density of 1.45g cm , which is 5.6 times higher than 
HATS-8b. This huge diversity in densities exhibited by 
super-Neptunes points to very different formation scenar¬ 
ios for these systems. In order to understand these differ¬ 
ences, however, it is incumbent on both precise ground- 
based surveys and future space-based transit surveys to 
expand the population of known super-Neptunes. 

We note that HATS-8b has a mass and period that 


puts it right o n the edg e of the ‘sub-Jupiter desert’ pro¬ 
posed bv lSzabo fc Kis s .2011;. This desert m ay be due 
to evaporation dKurokawa fc Nakamotol[2014f) . in which 
case HATS-8b is at the orbital limit for super-Neptunes, 
and if it were any closer to its host star then evaporation 
may have reduced it to a super-Earth class exoplanet. 
The fact that HATS-8 is currently a quiet star (see Sec¬ 
tion 12.41) may also point to a history of comparatively 
low X-ray flux, lowering the mean evaporation rate of 
HATS-8b and permitting it to exist on the edge of the 
“sub-Jupiter desert”. 


4.2. Ground-based Sensitivity 

While the period of HATS-8b is typical for an ex¬ 
oplanet detected by a wide-held, ground-based tran¬ 



Fig. 5.— The densities of exoplanets over four orders of mass. 
A characteristic “V” shape denotes two regimes of planets. The 
low mass planets (<~0.1Mj) show bulk densities that decrease 
with mass, as the gaseous planetary atmospheres grow in size rel¬ 
ative to the solid fraction of the planets. The high mass planets 
(>~0.1Mj) as gas planets where density increases with mass due 
to partial electron degeneracy in the planet interior. HATS-8b lies 
in the currently sparsely sampled transition region between these 
two regimes. All known exoplanets with well characterized densi¬ 
ties (density uncertainties less than 40%) are pl otted as blue circles , 
with data from the NASA Exoplanet Archive ijAkeson et al.|[2013D 
as of 2015 March 25. HATS-8b is plotted as red square. 

sit survey (P=3.583893 ± O.OOOOlOd), the derived mass 
(M=0.138± 0.019 Mj) and to a lesser extent the radius 
(R=0.873lg'j7 3 i?j), are much lower than typical, as is 
displayed in Figure[6] In fact HATS-8b is the third lowest 
mass planet to be detected by a wide-held, ground-based 
transit survey: o nly HAT-P-26b (IHartman et al.l 1201 lbT ) 
and HAT-P-llb (jBakos et al.l2010f are less massive. The 
discovery demonstrates that the HATSouth global net¬ 
work is capable of discovering gas planets in this low- 
mass regime. The importance of Keck/HIRES in this 
discovery should also be emphasized; given the magni¬ 
tude of HATS-8 and the mass of HATS-8b, Keck/HIRES 
is one of the only facilities in the world capable of con¬ 
firming this discovery. The ability to use a synthetic 
template (see Section IT41) for cross-correlation is also a 
significant development for faint targets such as HATS-8, 
for not only can it result in more precise radial velocity 
measurements, but it can also potentially dispense with 
the time-consuming step of creating an iodine-free tem¬ 
plate of the star. 

It is interesting to consider whether HATS-8b and 
other transiting exoplanets detected from ground-based 
transit surveys should be include d as targets on upcom¬ 
ing space missions such as T ESS (|Ricker et al.ll2014l) and 
PLATO (IRauer et al.l 120 ill) . Nominally HATS-8 would 
fall outside the magnitude range of the TESS mission 
(Ic < 12 for F,G,K stars), and so it would only be 
able to be monitored at a 30-minute cadence via the 
TESS full-framed images. However there are reason¬ 
able arguments to make that it would be wise to se¬ 
lect the 4 pixels which HATS-8 would occupy in order 
to obtain 2-minute cadence photometry. First, a set 
of well-sampled transits observed in 2019 would provide 
a long baseline with which to measure if the orbit had 
evolved in terms of transit timing variations. Secondly, 
the results from the Kepler and HARPS surveys for ex¬ 
oplanets show that many planetary systems are very 
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TABLE 5 

Stellar Parameters for HATS-8 


Parameter 

Value 

Source 

Identifying Information 

R.A. (h:m:s) 

19 h 39 m 46.08s 

2MASS 

Dec. (d:m:s) 

—25°44'53.9" 

2MASS 

2MASS ID 

2MASS 19394601-2544539 

2MASS 

Spectroscopic properties 

T eff * (K). 

5679 ± 50 

SPC a 

[Fe/H]. 

0.210 ±0.080 

SPC 

usin i (kms 1 ) ... 

2.00 ±0.50 

SPC 

7rv (kms- 1 ). 

l°g t? HK . 

19.958 ±0.041 

FEROS 

-5.18 ±0.1 

KECK 

Photometric properties 

v (mag). 

14.03 ±0.10 

NOMAD 

J (mag). 

13.098 ±0.024 

2MASS 

H (mag). 

12.779 ±0.029 

2MASS 

K a (mag). 

12.661 ±0.033 

2MASS 

Derived properties 

M* (M 0 ). 

1.056 ±0.037 

Isochrones±p*±SPC b 

R* (Ft©) . 

l.Uob_o .059 

Isochrones+p*+SPC 

log a * (cgs). 

4.386 ±0.071 

Isochrones+p* +SPC 

P* (cgs). 

1 l>5+ 0,21 

1 * 10 -0.35 

Isoc.hrones±p*±SPC c 

/.• (£©). 

1 11+0.31 

0.13 

Isochrones+p* +SPC 

My (mag). 

4.73 ±0.20 

Isochrones+p* +SPC 

Mk (mag,ESO) 

3.14 ±0.19 

Isochrones+p* +SPC 

Age (Gyr). 

5.1 ± 1.7 

Isochrones+p* +SPC 

Ay (mag) d . 

0.0000 ± 0.0094 

Isochrones+p* +SPC 

Distance (pc). 

829±“° 

Isochrones+p*+SP C 


a SPC = “Stellar Parameter Classification” method based on cross-correlating high- 
resolution spectra against synthetic templates ijBuchhave et al.l2012f) . These param¬ 
eters rely primarily on SPC, but have a small dependence also on the iterative anal¬ 
ysis incorporating the isochrone search and global modeling of the data, as described 
in the text. 

k Isochrones+p*+SPC = Based on the Y 2 isochrones \Yi et al.11200 if) , the stellar 
density used as a luminosity indicator, and the SPC results. 

c In the case of p* the parameter is primarily determined from the global fit to 
the light curves and radial velocity data. The value shown here also has a slight 
dependence on the stellar models and SPC parameters due to restricting the posterior 
distribution to combinations of p*+T e ff* + [Fe/H] that match to a YY stellar model. 
^ Total V band extinction to the star determined by comparing the catalog 
broad-band photometry listed in the table to the expected magnitudes from the 
Isochrones+p*+SPC model for the star. We use the lCardelli et al.| d 1989D extinction 
law. 


well aligned dFigueira et al.H2012t) , making it much more 
probable that additional planets could be discovered if 
we monitor a known transiting system. Although most 
hot Jupiters do no t appea r to have close planetary com¬ 
panions (jSteffen et al.ll20l3 ). it is possible that super- 
Neptunes such as HATS-8b would not follow this trend. 
Finally, precise photometry may help us understand this 
system in more detail, such as providing a stellar rota¬ 
tion rate or probing star spots if the planet happens to 
occu lt an active regi on of the star during transits (e.g. 
iMohler-Fischer et al.ll2013l) . 

4.3. Neighbors 

HATS-8 is in a relatively crowded field (see Table 0, 
that lies approximately 20° off the Galactic plane. While 
there is no evidence of a blended neighbor very close 
(< 5") to HATS-8, there are two neighboring stars at 
10" and 19.4", which are displayed in Figure [T] At the 
estimated distance to HATS-8 (8291+pc), these neigh¬ 
boring stars are not physically associated with HATS-8. 


They also do not significantly affect the HATSouth pho¬ 
tometry, as the HATSouth pixel scale (3.7" per pixel) 
means the neighbors are spatially resolved in the HAT¬ 
South images and in our follow-up photometry. How¬ 
ever the neighbor at 19.4" (2MASS19394689-2544386) 
is almost exactly the same magnitude (just 0.4 mag 
brighter) and with a similar color: (J — K) = 0.512 
compared with (J — K) = 0.437 for HATS-8. This makes 
2MASS19394689-2544386 an ideal reference star for long- 
slit, high resolution transit spectroscopy of HATS-8. 
Such observations have been carried out from ground- 
based facilities in other systems where t here is a n earb y 
reference star available (e.g. XO-2b. ISing et all 12012ft . 
We also expect HATS-8b to have an atmospheric scale 
height of 927 km if we assume a Neptune like mean molec¬ 
ular weight (^=2.53g.mol -1 ). Such a large scale height 
would increase the strength of the absorption signal for 
transmission spectroscopy. The primary limitation will, 
of course, be the fact that HATS-8 is relatively faint at 
V=14.03, so building up high signal-to-noise data will be 
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TABLE 6 

Parameters for the transiting planet HATS- 8 b. 


Parameter Value a 


Light curve parameters 

P (days) . 

T c (BJD) b . 

T 14 (days) b . 

T 12 = T 34 (days) b .. 

a/R* . 

C A’. ° . 

Rp/R * . 

b 2 . 

b = a cos i / R* . 

i (deg) . 


3.583893 ± 0.000010 
2456672.1102 ±0.0012 
0.1263 ±0.0030 
0.0109 ±0.0023 

q 91 +0.54 
y ' 2i -1.04 

17.42 ±0.32 

0.0814 ± 0.0038 

0 131+ 0 - 192 
u.ior_o.096 

r+0.21 


0.36Z 

87.8 


+1.2 


Limb-darkening coefficients d 


ci,i (linear term) . 0.2936 

C 2 ,i (quadratic term) .. 0.3212 

ci,r . 0.3906 

C2,r . 0.3083 


RV parameters 

K (ms -1 ) . 

e e . 

RV jitter (ms -1 ) f 

Planetary parameters 

M p (Mj) . 

Pp (Rj) . 

C(Mp,R p ) s . 

Pp (gem -3 ) . 

logffp (cgs) . 

a (AU) . 

Teq (K) h . 

© 1 . 

( F ) ( 10 9 ergs -1 cm -2 ) 1 


17.7 ± 2.5 
< 0.376 
0.00 ±0.99 


0.138 ±0.019 

0 - 87312 ;“® 

-0.08 

0.259 ±0.091 

9 

z.000 — o, 13 7 

0.04667 ± 0.00055 
1324I3® 
0.0138 ±0.0026 


6.94 


+ 1.82 


a The adopted parameters assume a circular orbit. Based 
on the Bayesian evidence ratio we find that this model is 
strongly preferred over a model in which the eccentric¬ 
ity is allowed to vary in the fit. For each parameter we 
give the median value and 68.3% (ler) confidence intervals 
from the posterior distribution. 

k Reported times are in Barycentric Julian Date cal¬ 
culated directly from UTC, without correction for leap 
seconds. T c : Reference epoch of mid transit that min¬ 
imizes the correlation with the orbital period. T 14 : to¬ 
tal transit duration, time between first to last contact; 
T 12 = T 34 : ingress/egress time, time between first and 
second, or third and fourth contact. 

c Reciprocal of the half duration of the transit used 
as a jump parameter in our MCMC analysis in place of 
a/R*. It is related to a/R* by the expression C,/R* = 
a/i?*(27r(l + e sin lj))/(P>/1 — b 2 y /1 — e 2 ) l|Bakos et al.l 

[2Ql0h . 

^ Values for a quadratic law, adopted from the tabu¬ 
lations by I Clare a 12004] ') according to the spectroscopic 
(SPC) parameters listed in Tabic [FI 

e The 95% confidence upper-limit on the eccentricity 
from a model in which the eccentricity is allowed to vary 
in the fit. 

^ Error term, either astrophysical or instrumental in ori¬ 
gin, added in quadrature to the formal radial velocity 
errors. This term is varied in the fit assuming a prior 
inversely proportional to the jitter. 

® Correlation coefficient between the planetary mass M p 
and radius R p determined from the parameter posterior 
distribution via C(M P , R p ) =< ( M p — < M p >)(R P — < 
R p >) > / (cMpCRp) > where < • > is the expectation 
value operator, and a x is the standard deviation of pa¬ 
rameter x. 

“ Planet equilibrium temperature averaged over the or¬ 
bit, calculated assuming a Bond albedo of zero, and that 
flux is reradiated from the full planet surface. 

1 The Safronov number is given by © = \ (Vesc /Vorb ) 2 — 
(a/R P )(M P /M ic ) (see lHansen Sz Barmanlf2007D . 

J Incoming flux per unit surface area, averaged over the 
orbit. 
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